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Using Lewis-Riesenfeld theory, we derive an exact non-adiabatic master equation describing the 
time evolution of the QED Schwinger pair-production rate for a general time-varying electric field. 
This equation can be written equivalently as a first-order matrix equation, as a Vlasov type integral 
equation, or as a third-order differential equation. In the last version it relates to the Korteweg-de 
Vries equation, which allows us to construct an exact solution using the well-known one-soliton 
solution to that equation. The case of time-like delta function pulse fields is also shortly considered. 
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I. INTRODUCTION 

Vacuum pair production by a strong electric field, pre- 
dicted by Schwinger in 1951 [ij, may now finally be seen 
due to the construction of ultra-strong laser systems 0. 
However, the corresponding fields are very different from 
the few special configurations for which an exact calcula- 
tion of the pair creation rate is possible. Thus, recently 
there has been increased interest in the development of 
approximation schemes, such as semiclassical methods 
[Jl^ and Monte Carlo simulations [Bj. 

A case that is relatively amenable to an exact treat- 
ment is the one of a purely time-dependent electric field. 
Here the spatial momentum is a good quantum num- 
ber, which allows one to reduce the time evolution of 
the system to a collection of mode equations labeled by 
the fixed momentum k. The pair production calculation 
can then be further reduced to a one-dimensional scat- 
tering problem, suitable for standard numerical or WKB 
methods 04^- Alternatively, the mode equation can be 
transformed to the quantum Vlasov equation, an integral 
equation for Afk{t), the total expected number of created 
pairs in the mode k [Tol - [l^ (see Ref. [lH| for a compari- 
son of the two approaches) . 

In this paper, we reconsider the time evolution of the 
QED Hamiltonian in a time-varying field using Lewis- 
Riesenfeld invariant theory [l6j and a suitable operator 
basis forming a spectrum generating algebra SU {1,1). 
We derive an exact non-adiabatic master equation for 
the time evolution of the Schwinger pair-production rate. 
This equation can be written equivalently as a first-order 
matrix equation, as a quantum Vlasov equation, or as a 
third-order differential equation. For a specific solution 
ansatz this third-order equation relates to the Kortcweg- 
de Vries (KdV) equation, which allows us to construct an 
exact solution using the well-known one-soliton solution 
to that equation. We also consider the case of alternating 
time- like delta function pulse fields, a type of fields which 



is of relevance for a recent proposal to apply Ramsey 
interferometry to the Schwinger effect [l7j . 



II. DERIVATION OF THE MASTER EQUATION 

We will give the derivation of the master equation for 
the scalar QED case; the derivation for the spinor QED 
case is similar, and will be included in a forthcoming, 
more detailed publication [l^. A scalar particle with 
charge q and mass m in a homogeneous time-dependent 
electric field with the gauge potential A\^{t) has the 
Fourier decomposed Hamiltonian of time-dependent os- 
cillators [in units of ft = c = 1] 



Hit) 



d?k 



(27r)3 



where 



uolit) = (fc|| -gyl||(t))2+kl+: 



(1) 



(2) 



We will quantize the theory in the Schrodinger picture, 
where the time-dependent quantum state obeys the func- 
tional Schrodinger equation 



(3) 



In this picture the field operators (^{x) and "kix) = t/)^ are 
time-independent with the momentum space commuta- 
tion relations 



(4) 



but the corresponding creation and annihilation opera- 
tors with the equal-time commutators 



(5) 
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are generally time-dependent |19l - l22| . 

TTk = dl{t)(pl{t) +b^k{t)(pk{t)- 



(6) 



Here ipk is an auxiliary field satisfying the classical mode 
equation 

yikit)+ujlit)^k{t)^0, (7) 

as well as the Wronskian constraint 

Wr[93fe, .^fe] = ifk{t)^*k{t) - ifl{t)ifk{t) = i. (8) 

The Eq. (O and the Wronskian determine Lpk{t) up to a 
phase factor, which we fix by requiring that Lpk{t) be real 
at the initial time ^o- Thus if io is finite, then for t < to 
one has 



,-iLJfc(0)(i-to) 



(9) 



(for to = —oo this should hold in the asymptotic sense). 
We note that the operators ak{t), and their hermi- 

tian conjugates are Lewis-Riesenfeld invariants, that is, 
they fulfill the Liouville-von Neumann equation 



i^^ + [ikit),Hk{t)]^0, 



(10) 



as can be easily checked. 

The ground state |Ofe,t) for the fc-th mode is annihi- 
lated by both ak{t) and b^k{t) and the n-th excited state 
is 



\nk,t} := j \Ok,t) 



nk'- 



(11) 



Thus the total time-dependent vacuum state is given by 
\0,t)=l[\Ok,t). (12) 

k 

In the free theory, the time-dependent vacuum state re- 
duces to the Minkowski vacuum, as expected. The scalar 
product for the quantized fields and their hermitian con- 
jugates allows us to find the Bogoliubov transformation 
between the past time to and the present time t, which 
is given by 

afe(^o) = f^k{to,t)ak{t) + i^k{to,t)b''_,,{t), 
blkito) - ^ll{to,t)bl,{t) + u*,{to,t)ak{t), (13) 



where 



Hk{to,t) = i'WT[(p*k{to),(pk{t)], 



(14) 



The Bogoliubov coefficients satisfy the relation for bosons 
\fJ.k{to,t)\'^ — Wk{to,t)\'^ = 1 . Our main object of interest. 



the mean number of pairs present at time t assuming 
that this number was Uk at the initial time can now 
be read off from 



{nk,t\al{to)ak{to)\nk,t) = \iyk{to,t)f{2nk + 1) + rik- 



Thus 



^k{t) := |i^fc(to,i)|'(2nfe + l) 



(15) 



(16) 



is the number of pairs spontaneously produced from the 
initial vacuum by the electric field. 

To obtain a time evolution equation for this quantity, 
we observe that the time-dependent Hamiltonian ([T]) has 
the spectrum generating algebra SU{1, 1). Choosing the 
Hermitian basis 



-Mf (to) - 



1 



(27r)3 
1 

(2^ 

i 

(2^ 



'al{to)ak{to) + b^kitoW-kito)], 
'ak{to)b-k{to) + al{to)blf^{to)], 
ak{to)b^k{to) - ai(to)^^fe(to)] , 



this algebra becomes 



(17) 



-2tM'^^\to). (18) 



The correlators are the expectation values of Eq. p7)) 
with respect to \nk, t) , that is, of the number of produced 
pairs and of pair creation and annihilation: 

l + 2Nk{t) = (2|!.fe(to,t)P + l)(2nfe + l), 



Mt\t) 



(pk(tQ,t)iyk{tQ,t) + iil{ta,t)vl{to,t)) 
x(2nfe + l), 

i{^ik{to,t)vk{to,t) ~ fil{tQ,t)vl{to,t)) 
x(2nfc + l). (19) 



Note that all three correlators are real and proportional 
to the quantum number 2nk + 1, and thus proportional 
to the ones defined by the vacuum state. 

Using Eq. (fT4|) and the mode equation ([7]) , we find the 
first order master equation 



a + 2Afk 



where 



/ dr^ 



n[r^ 
V -ni^^ 




l^kito) 



(20) 



(21) 
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with the initial conditions Nk — rikT M^j^^ = at t = 
(where may be — c»). An immediate consequence 
of the master equation (pOj) is the conservation of the 
quantity 

(1 + 2Mk? - {M^^^f - {M[7^f = (1 + 2n,)2. (22) 

This relates to the conservation of charge, as well as to 
the invariance of the Casimir operator for the SU{1, 1) 
algebra Ml- The spinor QED case can be treated anal- 
ogously |l8j ]. As far as the master formula (|20p is con- 
cerned, the generalization to the fermionic case requires 
only changing 1 -I- 2Afk{t) to 1 — 2J\fk{t), and replacing 
ivlit) by 



III. ALTERNATIVE FORMULATIONS OF THE 
MASTER EQUATION 

The first order matrix equation (j20p can be equiva- 
lently rewritten both as a single integral equation and as 
a third order linear differential equation. Since we work 
with a fixed mode k, in this section we will generally sup- 
press the index k and abbreviate loq := uj{t()). We will 
now also set Uk = 0. 

First, we combine the equations for Ai^^^ to a second 
order inhomogeneous equation for A^*^^\ 



dt 

dt 



+ (r!W)2x(-) 



(1 + 2J\f) 



(24) 



The homogeneous part of Eq. ((24| has the exact solutions 
X(-)(i)=C±e±^^4^*'"'^'(*') (25) 

with integration constants C^. Using those in the usual 
way to construct the solution of the inhomogeneous equa- 
tion with the appropriate initial conditions, we obtain the 
quantum Vlasov equation as the integral equation 



^{i + 2Af{t)) = n'^-\t) t dt' 

dt 



n'^-\t'){i + 2N{t')) 

(26) 



xcos(y di"f^(+)(i")) 



Second, inspection of the master equation ((20|) shows, 
that its general solution can be parameterized by a func- 
tion f[t) fulfilling the integral equation 



fit) 



Vl^~\t) 



Wo 



dt'f{t'){i^\t)+ic\t')) (27) 



with the initial condition /(to) = /(^o) = 0. Then, the 
correlators are given by 



l + 2Af = l+uJa dt'f{t')n^-\t'), 

J to 

X(-) = wofit), 

M^+^ = -coo f dt'f{t')n^+\t'). 

J to 



(28) 



Alternatively the integral equation ([27|) can, taking one 
derivative, be converted into a third order linear differ- 
ential equation. 



ujl{t) = {k\\-qA\\{t)f+tqEit)+kl+m^. (23) ^here 



F + Auj^F + 2{uj^yF - — ^ 



F{t) := f dt'f{t') 

Jtn 



(29) 



(30) 



and the initial conditions are i^(to) — F{to) — F{to) — 0. 
Observe that F is absent in Eq. p9p . which by Abel's 
theorem implies that the Wronskian of the solutions of 
the corresponding homogeneous equation is constant. 

The differential equation (I^Hl) bears an interesting re- 
lationship to the KdV equation. The form of the integral 
equation (P7)) suggests the ansatz 

/„)=<^, fw = ^!li>^. (31) 



8^0 



Defining r(t) := a;^(i)/a;o and then it(a;, i) :— —r{x—lOt), 
one can show that u satisfies the KdV equation. 



6uUx + Ut = 0. 



(32) 



Thus we can use certain solutions of the KdV equation to 
calculate pair creation rates for the corresponding electric 
fields. 



IV. EXACTLY SOLVABLE CASES 

We will now study two exactly solvable cases. First, we 
consider the following soliton-type solution of the KdV 
equation (see, e.g., Refs. |23l - l25| ) 



u{x, t) = —1 

which corresponds to 
2 



cosh^(x- lot) 



r{t) = 1 + 



cosh (ujot) ' 



4cjq cosh^(wot) 



(33) 



(34) 



This is a solution to Eq. ([29]) with the appropriate 
boundary conditions at to = — oo. The gauge potential 
is 



qA{t) = fcji - Jfcjf 



2w^ 



cosh (ujot) 



(35) 
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From Eq. ([28)1 we get the exact pair creation rate, 

1 



Scosh (woi) 



Note that Af{t) returns to zero for t — )■ oo, which is 
due to the sohtonic character that makes the scattering 
reflectionless. In fact, the mode sohition to Eq. ([7]) is 
given by 



-iojQt 



.A{t), 



each describing the exact time evolution of the cumula- 
tive pair creation variable A4 (t) for an electric field that 

(36) depends only on time, but is arbitrary otherwise. To the 
best of our knowledge, these equations are new. We have 
concentrated here on scalar QED, leaving the details of 
the spinor QED case to a more extensive publication [l8| . 

In future work, we also plan to study the precise con- 
ditions under which a non-adiabatic treatment is really 
necessary. To define the adiabatic approximation, we 
write the mode solution in terms of the adiabatic basis 

(37) M 



where the amplitude is 

A{t) = (e^"^"* + l)^2i^i(2, 2-i;l-i; -e^""*) (38) 

with 2F1 the hypergeometric function, and it does not 
have a negative frequency part in the future. The Bo- 
goliubov coefficient ([Ti)) is 



Ki) 



2^0 



-A*{t), 



(39) 



which approximately leads to JV{t) = 2e^""* for uqI <C 
— 1 and Af{t) = 26"^""* for ujot ^ 1, the leading approx- 
imation to the exact formula (j36p . Thus there is no pair 
creation in this case, contrary to the somewhat similarly- 
looking Sauter field case [2a|. This example also shows 
clearly that, as emphasized in Ref. no direct physi- 
cal meaning should be ascribed to A4 (t) at intermediate 
times. 

Second, we consider an electric field consisting of two 
opposite delta function pulses. 



(40) 



which has the gauge potential of a potential well [17| : 
A|| = for t < and t > ti, corresponding to u}k{0), 
and A|| — —Eq for < i < ii, corresponding to uJk- The 
master equation (l20t together with continuity at t = 
leads to the pair production for the period < t < ti 



l + 2A4(i) = (2nfe + l)(^ 



1 - 



2ujk 



(+) 



cos(2a;fci) , (41) 



and for the period t > ti it now remains constant, retain- 
ing its value for ti. Note that for a single delta function 
pulse Mk{t) keeps oscillating, so that the limit t — ?> 00 
cannot be defined. This is presumably due a combina- 
tion of the unphysical character of such a field and the 
non-Markovian nature of the time evolution. 



-ie(t) 



-Mt) 



(42) 



where 9{t) = J^^ dt'ujk{t') and the Bogoliubov relation 
IttfcP — IPkl"^ = 1 holds. We then replace \i^k{to,t)\'^ by 
|/3fc(i)P in the definition (HH) oiJ\f{t). 

From Eq. (fT4|) one can easily show that for this ap- 
proximation to hold it is sufficient to assume that 



Wfc(O) V '^kit) 





LOkit) 


1 





« mt)\ (43) 



throughout the time evolution. This criterium is similar, 
although not strictly equivalent, to the one given in [ll| . 



< 1, 



< 1. 



(44) 



In any case, all the inequalities in (j43| . (|44l) are certainly 
fulfilled for even the strongest laser sources which are 
presently existing or in development. Those have a field 
strength still much lower than the critical strength Ec = 
m^/e and the characteristic time scale much longer than 
the Compton time 

Concerning the relation of the master equation to the 
KdV equation, although there is a well-known connection 
between the latter equation and one-dimensional quan- 
tum mechanical scattering (see, e.g., Refs. [23l - l25i |27| ) . 
it appears not to have been previously applied to the 
Schwinger pair creation problem. It will be interesting to 
see whether also the multi-soliton solutions of the KdV 
equation may be used in this context. 

Finally, let us mention that it is straightforward to 
extend our master equation to the case of an initial state 
which is a thermal state at temperature T. As will be 
shown in Ref. [31, such a change leads again only to an 
overall factor (coth(/3wfe(0)/2) -I- l) multiplying all three 

correlators 1 -I- 2A4, so that the master equation 

itself remains unaffected. 
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